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Using form factors enrolled to the transition matrix elements and calculated via 
light-cone QCD sum rules including next-to-leading order corrections in the strong 
coupling constant, we provide a systematic analysis of the Bg — )• fo{980)£'^ £~ both 
in the standard and universal extra dimension models. In particular, we discuss 
sensitivity of the differential branching ratio and various double lepton polarization 
asymmetries on the compactification factor of extra dimension and show how the 
results of the extra dimension model deviate from the standard model predictions. 
The order of branching ratio makes this decay mode possible to be checked at LHCb 
in near future. 
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I. INTRODUCTION 



It is well known that the decays of the Bg meson are very promising tools to constrain 
the standard model (SM) parameters, serve to determine the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, enable us to understand the origin of the CP violation 
and help us search for new physics (NP) effects beyond the SM. One of the possible decay 
modes of the Bg meson is the semileptonic Bg — )■ fo{980)i~^i~ . This channel proceeds via 
flavor- changing neutral currents (FCNC) transition of 6 — )■ si^i^ at quark level, which 
is induced at loop level in the SM and is therefore sensitive to NP effects. The extra 
dimensions as NP effects can contribute to such loop level transitions and enhance the 
branching ratio. It has been shown that, in the presence of a single universal extra dimension 
(UED) compactified on a circle with radius R, the branching ratios of the Bg — )■ r]^'H~^ i~ / uu 
and Bg (puv which are also based on 6 — )• st^l~ jvv^ increase significantly at lower values 
of the compactification factor, [ll, 

In the present work, we investigate the effect of the UED on some physical observables 
related to the semileptonic Bg — )• /o(980)£"'"£~. The UED with a single extra dimension 
called the Appelquist, Cheng and Dobrescu (ACD) model is a kind of extra dimension 
(ED) 0^ which allows the SM fields (both gauge bosons and fermions) to propagate in the 
extra dimensions (for more details see for instance The ACD model has been previously 
applied to many decay channels. For some of them see [8-2l| and references therein. 



In the SM, the effective Hamiltonian describing the Bg — )■ /o(980)£"'"£ transition at quark 
level can be written as: 



GFOiemVtbVf*g 



(1) 



where Vij are elements of the CKM matrix, aem is the fine structure constant, Gi? is the Fermi 
constant, and Cj^^ , Cg-^^ and Cio are Wilson coefficients. In the ACD model, the form of 
the effective Hamiltonian remains unchanged but due to the interaction of the Kaluza-Klein 
(KK) particles with the usual SM particles and also with themselves, the Wilson coefficients 



are modified. This modification is done in [22|-|26| in leading logarithmic approximation in 



a way that each Wilson coefficient is written in terms of some periodic functions as: 

00 

F{xt, 1/R) = Fo{Xt) + Fn{Xt, Xn), (2) 

n=l 

where, -Fo(xj) is ordinary SM part and the other part can be written in terms of the com- 
pactification factor 1/R, by means of the following definitions: 

Xt = m'^jM^, Xn = ml/m'^^, mn = n/R, (3) 

where, rrit is the mass of the top quark, Mw is the mass of the W boson, rrin is the mass 
of the KK particles and n = corresponds to the ordinary SM particles. Few comments 
about the lower bound of the compactification factor are in order. From the electroweak 
precision tests, the lower limit for 1/R had been previously obtained as 250 GeV in js], 
if Mh > 250 GeV expressing larger KK contributions to the low energy FCNC processes, 
and 300 GeV if Mh < 250 GeV. Analysis of the B — > ^^7 transition and also anomalous 
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magnetic moment had shown also that the experimental data are in a good agreement 
with the ACD model if 1/R > 300 GeV 23]. Taking into account the leading order (LO) 
contributions due to the exchange of KK modes as well as the available next-to-next-to- 
leading order (NNLO) corrections to also B{B — )• Xg'j) in the SM, the authors of [2l| 
have obtained a lower bound on the inverse compactification radius 600 GeV. Using the 
electroweak precision measurements and also some cosmological constraints, the authors of 



28[ and [29[ have found that the lower limit on compactification factor is in or above the 
500 GeV range. We will plot the physical observables under consideration in the range 
1/R E [200 — lOOOjGeT^ just to clearly show how the results of the UED deviate from those 
of the SM and grow decreasing the 1/R . 

The numerical values for different Wilson coefficients both in the SM and ACD model in 
the range 1/R E [200 — lOOOjGel^ are presented in Table [H From this Table, we see that 
in the ACD model, the Giq is enhanced and G^^^ is suppressed considerably in comparison 
with their values in the SM. 



1/R [GeV] 




Cio 


^9 


200 


-0.195212 


-5.61658 


4.83239 + 3.59874i 


400 


-0.266419 


-4.65118 


4.7538 + 3.54366i 


600 


-0.283593 


-4.43995 


4.7366 + 3.53161i 


800 


-0.29003 


-4.36279 


4.73032 + 3.52721i 


1000 


-0.293092 


-4.32646 


4.72736 + 3.52514i 


SM 


-0.298672 


-4.26087 


4.72202 + 3. 52139i 



TABLE I. Numerical values for Cj , Ciq and values of Cg at transferred momentum square, 
q'^ = 14 for different values of 1/R as well as the SM 



Here we should mention that, besides the aforementioned contributions, the Wilson co- 
efficient Gg^^ receives also long distance contributions from J/ ip family parameterized using 



Breit-Wigner ansatz [30 



I.e., 



3n _ ^ TiV-^i^n^v, 



where, G^^^ = 0.362. As the phenomenological factors, sei have not known for the transition 
under consideration, we will choose the values of the which do not lie on the J/ ip family 
resonances. This is possible in the case of the differential branching ratio and double-lepton 
polarization asymmetries under consideration in this work, however, to calculate the total 
branching ratio as well as the average double-lepton polarization asymmetries, which require 
integration over g^, one should take also into account such contributions. For more details 
about the long distance contributions see for instance j3l|. 



The layout of the paper is as follows. In the next section, we present the transition 
matrix elements expressed in terms of form factors and the formula for decay rate as well as 
various double lepton polarization asymmetries. In the last section, we numerically analyze 
the observables in terms of compactification factor, 1/i? of the ACD model and discuss the 
results. 
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TABLE II. Bg — )• /o(980) transition form factors including next-to- leading order corrections 33l |. 

F^{q^ = 0) a, bi 



Fi 0.238 ±0.036 l.50tom O-SStao? 
Fo 0.238 ± 0.036 0.53l|j:J^ -0.36t[{;[j^ 
Ft 0.308 ±0.049 1.461^;^^ 0-58lo:o7 



II. TRANSITION MATRIX ELEMENTS AND OBSERVABLES RELATED TO 

THE Bs /o(980)£+^' CHANNEL 

To find the amplitude of the decay channel in question, we need to sandwich the afore- 
mentioned effective Hamiltonian between the initial and final states. As a result, we obtain 
the following transition matrix elements defined in terms of form factors Fo(g^), Fi[q^) and 



(/ob/o)IW5&|i?.(PBj) = -^{Fl{q' 



P„ 



m 



— m 



fo 



+ Fo{q' 



mi 



m 



fo 



}, (4) 



{fo{.Pfo)\s(^i^yl^q''^Bs[pBs)) = ^— g^^M - i'^Bs - ' (5) 

where P = ps^ + Vfo q = pb^ ~ Pfo- -^^^ simplicity in some parts of calculations, it is 
convenient to introduce the auxiliary form factors /+ and /_, 

(/o(p/o)I^7m75&|5,(pbJ) = -^ {Uiq')P, + f-{q')q,} 

(6) 



such that 



Fi{q') = U{q') 



q 



2 



Fo{q') = IW) + 2 2 /-(g^) • (7) 



The form factors, Fi, Fq and Ft are calculated via light-cone QCD sum rules both at the 
leading order and the next-to- leading order corrections in [33|]. We use the latter to analyze 
the considered physical observables. The fit parametrization of the form factors including 



next-to-leading order corrections in as is given as [33 



Fm_ 

1 - aiq^/m\^ + bi{q'^/m%J' 



^iil^) - 1 „ „2/„„2 * , I. /„2/„,2 N2 ' (^) 



where Fi denotes any function among -Fi_o,t- The parameters Oj and hi as well as -Fj(O) are 
given in Table [Tll 

Now we proceed to calculate some observables such as differential decay rate and double 
lepton polarization asymmetries. With the matrix elements in terms of form factors one can 
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easily obtain the dependent differential decay rate as: 



-{q\l/R)- 



3g2 



Qmi\C,,{llR)\\m's^-miyFi{q') 



+ + 2ml)\ 



C^/\q^l/R)F,{q') + 



2C'j"il/R)im,-m,)FTiq' 



+ \C,o{l/R){'{q'-^mi)XF^{q'' 



(9) 



with f = Y 1 A = A(m^^, m^^j, g^) with A(a, b,c) = {a — b — c)^ — 46c and rrii is the 

lepton's mass. 

To calculate the double-polarization asymmetries, we consider the polarizations of both 
lepton and anti-lepton, simultaneously and introduce the following spin projection operators 
for the lepton i~ and the anti-lepton (see also 34-36|): 



Ai=^(i + 75/r) 

A2 = ^(l + 75/+) 



(10) 



where i = L,N and T correspond to the longitudinal, normal and transversal polarizations, 
respectively. Now, we define the following orthogonal vectors in the rest frame of lepton 
and anti-lepton: 



'AT 



0, Crp 



0,e + 



0,62^ 





0, 



P- 



' \p-\J 

Pfo X P~ 



\Pfo X P 

0, Cjy X 
P+ \ 



0, 
0, 



(0, 



\p+\J 

Pfo X P+ 

\Pfo X P+l 



X 



(11) 



where p^p are the three-momenta of the leptons £^ and is three- momentum of the final 
/o meson in the center of mass (CM) frame of i~ . By Lorenz transformations, the 
longitudinal unit vectors are boosted to the CM frame of l~t^, 
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while the other two vectors are kept the same. Now, we define the double-lepton polarization 
asymmetries as 34j-l36| : 







dr 




4) 


ds 


dT^^ 




dT 

ds 



'-Si ,s^) 



dV dV 



(13) 



'^i 1 Sj J 



where the subindex j also stands for the L, N or T polarization. The subindexses, i and 
j correspond to the lepton and anti-lepton, respectively. Using the above definitions, the 
various 1 / /^-dependent double lepton polarization asymmetries are obtained in the following 
way: 



3A(s, 1/R) 

,2 ^2 



-12mj,mfs\D\' + mj,JC\'i2X - (1 - v')i2X + 3(1 - f^J^))], 



sA{s,l/R) 



Pll{sMR) 

Pln{s, l/R) 
PndsMR) 

^-(^'^/^) sA^sMR) 
Ptl{sMR) = Plt{sMR) 

Pmt{s, l/R) 

PTNis, l/R) 

PtAs, l/R) 

PNNis, l/R) 



^= Re[-2Am\m1{l - ff^)C*D + A'm|^(l + v^)\A\ 



3A(s,l/i?) 
-PnAsA/R), 



Re[mB,miv{l - r/J|C|^ + mB,rhivsC*D], 



Im[2X m%^A*C] 



(14) 

(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



Am, 



Re[-2Amlm^{l - rf„)C*D - A'm|^(l + v^)\A\^ - l2mlm^s\D\ 



3A{s,l/R) 

+ m|JCp{2A' - (1 - v'){2X' + 3(1 - f^J^)}], 



(21) 



4m 



3A{s,l/R) 



Re[2Amlmf{l - rfy)C*D - A'm|^(3 - v^)\A\'^ + l2mlm^s\D\ 



+ mlJCn2X -{l-v'){2X -3{l-rf,f)}], 



(22) 



where, s 



mt, 'JO m- 



= ^' = ^(l'^/o>^) and 



A{s,l/R) 



4m 



Bs 



Re[2Am%mf{l - ffo)D*C + X m%^{3 - v^)\A\^ + l2m%mfs\D 



+ m|JCp{2A' - (1 - v'){2X' - 3(1 - rj^f)}], 



(23) 



with 



A = A{s, l/R) = 2Cl"{s, l/R)f+{s) - 4C'/^{l/R){mh + m. 



Fris) 



B = B{s, l/R) = 2Cl^^{s, l/R)f-{s) + AC'j"il/R)imb + m,) 

C = Cis,l/R) = 2C,oil/R)Uis), 
D = D{s,l/R) = 2Cio{l/R)f-{s) . 



effl 



^Bs + rrif^ ■ 
Ft{s) 



[rriBs +m/Jsm|^ 



m 



Bs 



m 



(24) 
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III. NUMERICAL RESULTS 

In this section, we numerically analyze the physical observables and discuss their sensitiv- 
ity to the compactification factor of extra dimension. The main input parameters are form 
factors in the matrix elements whose fit parametrization are presented in the previous sec- 
tion. To proceed in numerical calculations, we also need to know the numerical values of the 
other input parameters. We use the values: rrit = 167 GeV, mw = 80.4 GeV, nih = 4.8 GeV, 
rris = 0.14 GeV, = 0.105 GeV, = 1.778, \VtbV*\ = 0.041, Gp = 1-166 x 10"^ Gel/-^, 
"em = ily, = 1.42 X lO^^^ ^ ^ 0.980 GeV and rriB, = 5.36 GeV. 

Considering the central values of the form factors, we plot the dependence of the differen- 
tial branching ratio and various double lepton polarization asymmetries for the Bs — ?• foC-'^^~ 
decay channel on the compactification factor (1/-R) of the extra dimension in figures 1-8. As 
the results of the electron channel are very close to those of the /x, we will depict only the 
results of the and r channels. As it is evident from the formulas in the previous section 
that the observables depend on s, we will present our results at three fixed values of this 
parameter for the /x and two fixed values for r channel in the allowed kinematical region 
[Ami s < [1 — \J^oY)- Note that in each figure we see graphs of the lines with the same 
colors. The straight line in each case depicts the result of the SM and the curve line stands 
for the ACD model prediction. 




200 400 600 800 1000 

l/R[GeV] 



FIG. 1. Dependence of the branching ratio on the \/R for muon channel at three fixed values of 
the s. The blue, green and red lines belong to the values s = 0.2, s = 0.3 and s = 0.5, respectively. 
The straight line shows the result of SM and the curve depicts the ACD model prediction in each 
case. 

From these figures, we obtain the following results: 

• There are considerable discrepancies between the results of the UED and SM pre- 
dictions at lower values of the compactification factor for all observables and both 
lepton channels. When 1/i? is increased, the differences between the predictions of 
two models become small so that two models have approximately the same predictions 
at 1/i? = 1000 GeV. 

• As it is expected, the branching ratio in r channel is small comparing to that of the 

• An increase in the value of s ends up in a decrease in the value of the differential 
branching ratio. 
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200 400 600 800 1000 

l/R[GeV] 



FIG. 2. Dependence of the branching ratio on the \/R for tau channel at two fixed values of the 
s. The blue and green lines belong to the values s = 0.5 and s = 0.6, respectively. The straight 
line shows the result of SM and the curve depicts the ACD model prediction in each case. 



0.90 




^ ^ 0.55 C • • • ' • • • ' • • • ' • • • 3 

200 400 600 800 1000 200 400 600 800 1000 

l/R[GeV] l/R[GeV] 

FIG. 3. Dependence of the Pll and Pnn on the 1/R for muon channel at three fixed values of the 
s. The blue, green and red lines belong to the values s = 0.2, s = 0.3 and s = 0.5, respectively. 
The straight line shows the result of SM and the curve depicts the ACD model prediction in each 
case. 

• The deviations of the UED results from those of the SM on double lepton polarization 
asymmetries are small in comparison with the deviation of the differential branching 
ratios from corresponding SM values. However these can not be overlooked. 

• The contributions of KK modes enhance the absolute values of the Pll for r as well 
as the PnN) Ptt and Plt for both lepton channels, but they decrease the absolute 
values of the Pll for ix and Pln and Ptn for both leptons. 

• The Pll for ^ and Ptt for r have negative signs but the rest of double lepton polar- 
ization asymmetries have positive signs. 

Now, we would like to discuss how the uncertainties of the form factors affect the physical 
quantities under consideration. For this aim, we plot the aforementioned physical observables 
on the compactification factor in figures 9-16 when the uncertainties of the form factors are 
taken into account. These figures are plotted at s = 0.2 and s = 0.6 for the fi and r channels, 
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FIG. 4. Dependence of the Pll and Pnn on the 1/R for tau channel at two fixed values of the s. 
The green and red lines belong to the values s = 0.5 and s = 0.6, respectively. The straight line 
shows the result of SM and the curve depicts the ACD model prediction in each case. 



0.85 




0.60 



200 400 600 800 1000 200 400 600 800 1000 

l/R[GeV] l/R[GeV] 



FIG. 5. The same as FIG. 3 but for Ptt and Pln- 




400 600 800 

l/R[GeV] 



J 0.07 I — 

1000 200 



600 
l/R[GeV] 



FIG. 6. The same as FIG. 4 but for Ptt and Pln- 



respectively. From figures 9 and 10 for the branching fractions, it is clear that the difference 
between the UED and SM models predictions exist and can not be killed by uncertainties 
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FIG. 8. The same as FIG. 4 but for Plt and Ptn- 




200 400 600 800 1000 

l/R[GeV] 



FIG. 9. Dependence of the branching ratio on the \/R for muon channel at s = 0.2 when errors of 
the form factors are taken into account. The straight band shows result of the SM and the curve 
band refers to the ACD model prediction. 

of the form factors especially at lower values of the compactification factor. The figures 
11-16 for double lepton polarization asymmetries also depict that except the Ptt and Ptn 
in r channel and Plt in both lepton channels, there are discrepancies between two model 
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200 400 600 800 1000 

l/R[GeV] 



FIG. 10. Dependence of the branching ratio on the 1/R for tau channel at s = 0.6 when errors of 
the form factors are taken into account. The straight band shows result of the SM and the curve 
band refers to the ACD model prediction. 




200 400 600 800 1000 200 400 600 800 1000 

l/R[GeV] l/R[GeV] 

FIG. 11. Dependence of the Pll and Pnn on the 1/R for muon channel at s = 0.2 when errors of 
the form factors are taken into account. The straight bands show results of the SM and the curve 
bands refer to the the ACD model predictions. 

predictions at lower values of the 1/R even if the errors of the form factors are encountered. 

In conclusion, making use of the related form factors calculated via light-cone QCD 
sum rules up to next-to-leading order corrections in a^, we analyzed the sensitivity of the 
differential branching ratio and various double lepton polarization asymmetries on the com- 
pactification factor of extra dimension for the Bg — )■ /o(980)£"'"£~ transition. Our numerical 
calculations depict considerable deviations of the extra dimension model results from the 
SM predictions. These differences can not be killed by errors of the form factors in the 
allowed regions of the compactification parameter previously discussed. Such discrepancies 
can be interpreted as signals for existing extra dimensions in nature which can be searched 
for at hadron colliders. As a final note it is worth to estimate the accessibility to measure 
the branching ratio and lepton polarization asymmetries. An observation of a 3 cr signal for 
asymmetry of the order of the 1% needs about ~ 10^^ BB pairs. This allow us to measure 
the branching ratio and the polarization asymmetries shown in the figures (1)-(16) in princi- 
ple. The order of branching ratios show that the Bg — )■ /o(980)^"'"£~ decay channel both for 
11 and r leptons can be detected at LHC. However, as experimentalists say, there are some 
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FIG. 12. Dependence of the Pll and Pnn on the 1/R for tau channel at s = 0.6 when errors 
of the form factors are taken into account. The straight bands show results of SM and the curve 
bands refer to the ACD model predictions. 
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FIG. 13. The same as FIG. 11 but for Ptt and Pln- 
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FIG. 14. The same as FIG. 12 but for Ptt and Pln- 



technical difficulties to measure the lepton polarizations. In the case of jj, this lepton should 
be stopped in order to measure its polarizations which is not yet possible experimentally. 
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200 400 600 800 1000 200 400 600 800 1000 

l/R[GeV] l/R[GeV] 



FIG. 15. The same as FIG. 11 but for Plt and Ptn- 




200 400 600 800 1000 200 400 600 800 1000 

l/R[GeV] l/R[GeV] 



FIG. 16. The same as FIG. 12 but for Plt and Ptn- 

For r lepton, we should reconstruct then analyze the decay products of this lepton. In this 
case, we face with the problem of the efficiency of the reconstruction. If these technical 
difficulties over come, by measuring the considered double-lepton polarization asymmetries, 
we can get valuable information about the nature of interactions included in the effective 
Hamiltonian because as the large parts of the uncertainties are canceled out, the ratio of 
physical observables such as CP, forward-backward asymmetry and single or double-lepton 
polarization asymmetries less suffer from the uncertainty of the form factors compared to 
the branching ratio. 
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